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Determinations of a s using JADE data of e + e Annihilations 
at y/s = 22 to 44 GeV 

P.A. Movilla Fernandez a 

a III. Physikalisches Institut A, RWTH Aachen — Physikzentrum, 
Sommcrfcldstr., D-53056 Aachen, Germany 

Data recorded by the JADE experiment at the PETRA e + e~ collider were used to measure distributions of 
new event shape observables. The distributions were compared with resummed QCD calulations (C(Qg)+NLLA), 
and the strong coupling constant a s was determined at y/s — 22, 35 and 44 GeV. The results are in agreement 
with previous combined results of PETRA but have smaller uncertainties. Together with corresponding data from 
LEP, the energy dependence of a s is significantly tested and is found to be in good agreement with the QCD 
expectation. 



1. INTRODUCTION 

Summaries of measurements of a s from various 
processes and at different energy scales Q demon- 
strate that the energy dependence of a s is in 
good agreement with the prediction of Quantum 
Chromodynamics (QCD). High precision tests are 
provided by e + e _ annihilation data which were 
collected between some GeV and nearly 200 GeV. 
Unfortunately, in the interesting intermediate en- 
ergy region where a s runs very quickly, the "best" 
value of a s , 0.14 ± 0.2 §, is a combined result at 
y/s = 35 GeV based on measurements at PETRA 
that are affected by large theoretical uncertainties 
because only QCD predictions up to 0{p? s ) were 
available for event shape observables. 

Significant progress has been made in pertur- 
bative QCD calculations since 1992. Observables 
have been proposed for which perturbative pre- 
dictions are extended beyond the next-to-leading- 
order through the inclusion of leading and next- 
to-leading logarithms || which are summed to all 
orders of a s (NLLA). Further, the precise mea- 
surements of the hadronisation process have im- 
proved the modelling of hadronic final states with 
Monte Carlo programs. 

It is desirable to apply these LEP-established 
developments consistently to lower energy data 
in order to improve quantitative studies of the 
running of a s (Q) over large ranges of the energy 
scale Q, using identical experimental techniques 



and theoretical calculations in order to minimise 
point-to-point systematic uncertainties. In the 
following, ©(cv^+NLLA determinations of a s us- 
ing data from the JADE experiment at PETRA 
at y/s — 22, 35, and 44 GeV are presented 

2. JADE DATA 

The JADE detector @ was one of the five ex- 
periments at the PETRA electron-positron accel- 
erator. It has operated from 1979 until 1986 at 
centre-of-mass energies of y/s = 12 to 46.7 GeV. 
The main components of the detector were the 
central jet chamber to measure charged particle 
tracks and the lead glass calorimeter to measure 
electromagnetic showers. 

Multihadronic events were selected by the stan- 
dard JADE selection cuts fl which reduce the 
background from 77 and r-pair events to less than 
0.1% and 1%, respectively (i). The final numbers 
of events which were retained for this analysis are 
listed in Table [j], together with the correspond- 
ing integrated luminosities. The physical results 
of previous publications by JADE can be very 
well reproduced 

Monte Carlo detector simulation correspond- 
ing to the JADE data were retrieved for 35 and 
44 GeV. They were generated using the JETSET 
QCD shower event generator version 6.3 fl. In 
general, there is good agreement of the detector 
simulation with data for all event shape distribu- 
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Table 1 

Number of multihadronic events in JADE data and 
corresponding integrated luminosities. 



year 


[GcV] 


data 


luminosity 


1981 


22 


1404 


2.4 pb- 1 


1984/85 


39.5-47.8 


6158 


80 pb- 1 


1986 


35 


20 926 


40 pb- 1 



tions studied here. The simulated data can thus 
be used to correct for detector effects in the mea- 
sured data. 

3. DETERMINATION OF a s 

3.1. Event shapes and jetrates 

From the data samples, the event shape distri- 
butions of thrust T, the heavy jet mass Mjj, the 
total and wide jet broadening Bt and Byy and 
the differential 2-jet event rate using the Durham 
jet finder j/23 were determined (c.f. g). For the 
latter three observables, no experimental results 
have previously been presented at these energies. 

3.2. Correction procedure 

The event shape data were corrected for the 
limited acceptance and resolution of the detector 
and for initial state photon radiation effects by 
applying a bin-by-bin correction procedure. It 
turns out that there is an excellent agreement 
between the data and the model prediction at 
hadron level over the whole kinematic range of 
the observables even for those ones which were 
never tested at such low energies Qj. 

No detector simulation were available for the 
22 GeV data. Therefore only the differential 2- 
jet rate was considered because it is known to de- 
pend to a lesser extend on detector effects. The 
35 GeV correction was applied to the 22 GeV 
data since the correction factors for this observ- 
able were found to show no significant energy de- 
pendence [||. 

In a second step, the data distributions were 
corrected for hadronisation effects by applying 
bin-by-bin correction factors derived from the ra- 
tio of the JETSET generated distributions before 
and after hadronisation. The data distributions, 
thus corrected to the parton level, can be com- 
pared directly to analytic QCD calculations. 
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Figure 1. Distributions for 1 — T and Mh / '\/s at 
-y/s = 35 GeV corrected to parton level. The fits of 
the QCD predictions are overlaid and the fit ranges 
are indicated by the arrows. 

3.3. 0(a^)+NLLA predictions 

The event shape observables used in the analy- 
sis are predicted in perturbative QCD by a combi- 
nation of the 0(a 2 s ) @] and the NLLA g calcula- 
tions. The 0(a 2 s ) calculation yields an expression 
of the form 

R 0al {y) = l + A{y)[^j+B{y)(^j , 

where R(y) = Jq dy 1/co • da/dy is the cumula- 
tive cross-section of an event shape observable y 
normalised to the lowest order Born cross-section 
(To . The NLLA calculations give an expression for 
R(y) in the form: 

W(y) = (x + c 1 (|) + c 2 (g) 2 ). 
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Figure 2. Distributions for Bt and Bw at y/s = 
35 GeV corrected to parton level. The fits of the 
QCD predictions are overlaid and the fit ranges are 
indicated by the arrows. 
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Figure 3. Distributions for 2/23 at \/s = 22 and 
35 GeV corrected to parton level. The fits of the 
QCD predictions are overlaid and the fit ranges are 
indicated by the arrows. 
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where L — ln(l/y). The functions g\ and gi are 
given by the NLLA calculations. The coefficients 
C\ and Ci are known from the 0{a 2 s ) matrix el- 
ements. 

The strong coupling a s was determined by x 2 
fits to the event shape distributions corrected to 
the parton level. For the sake of direct com- 
parison with other published results this analysis 
closely followed the procedures described in ||. 
The so-called ln(i?)-matching scheme was chosed 
to merge the 0(a 2 s ) with the NLLA calculations. 
The renormalisation scale factor, x^ = /u/\/s, was 
set to x^ — 1 for the main result. 

The fit ranges are determined by the require- 
ments of moderate hadronisation uncertainties 



within the fit ranges, good x 2 /d.o.f. and reason- 
ably stable fits under changes of the range. 

3.4. Systematics 

To study systematic uncertainties details of 
the event selection and of the correction proce- 
dure were modified. The impact of the hadroni- 
sation model of the JETSET generator was stud- 
ied by varying several significant model parame- 
ters. Mass effects are estimated conservatively by 
excluding the b-quark from the correction proce- 
dure. Uncertainties due to contributions of un- 
known higher order corrections are estimated by 
the changes of a s when varying the renormali- 
sation scale factor x^ in the range of 0.5 to 2.0. 
Any resulting deviation from the main result was 
considered as systematic error. 
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3.5. Results 

As an example, the theoretical predictions for 
the observables 1 — T, Mjf, and at 
yfs = 35 and ?/23 at yfs = 22 and 35 GeV are pre- 
sented in Figures [j] - 1| and superimposed on the 
corresponding measured distributions corrected 
to partem level. In general, resummed calcula- 
tions (solid lines) give a good description of the 
observables over their entire range of values even 
at the low PETRA energies, while the G{a 2 s ) cal- 
culations (dashed lines) which are shown for com- 
parison fail in describing the 2-jet region. 

The values of a s and the errors obtained at 35 
and 44 GeV are shown in Figure || in compar- 
ison with the a s values measured by OPAL at 
y/s = M z o JH) . The values of a s exhibit a similar 
scattering pattern at all energies. This demon- 
strates the strong correlation of the systematic 
uncertainties at JADE and OPAL. 

The individual results of the five observables at 
each energy were combined into a single value. 
The final results for a, are 



0.0069 



(syst.) 



a s (44GeV) = 0.1372 ± 0.0017(stat.)I 
a s (35GeV) = 0.1434 ± 0.0010(stat.)+°;°Jg5(syst.) 
a s (22GcV) = 0.1608 ± 0.0083(stat.)+o °o5g(syst.) 

where the result at 22 GeV is based on y 2 3 only. 
The main contributions to the systematic er- 
rors come from renormalisation scale uncertain- 
ties and mass effects which causes systematically 
higher values for a s . 

4. SUMMARY AND CONCLUSIONS 

JADE data at y/s = 22, 35 and 44 GeV were 
used to measure a s simultaneously with old (T, 
Mh) and new (Bt, %, 2/23) event shape observ- 
ables. The corresponding 0(a;?)+NLLA predic- 
tions are in good agreement with the measured 
distributions corrected to parton level even at 
PETRA energies. 

It can be pointed out that the present measure- 
ments of a s are the most precise in this energy 
region of the e + e _ continuum. The individual re- 
sults are consistent with each other. The errors 
are dominated by effects of missing higher orders 
in the QCD calculations. 



OPAL 91 GeV 



JADE 35 GeV 



T 



D, 



weighted average h -"*-"-> — -' 
0.08 0.1 0.12 



0.14 



0.16 



a. 



OPAL 91 GeV JADE 44 GeV 



T 



D, 



weighted average h 



0.08 



0.1 



0.12 



0.14 



0.16 



0t„ 



Figure 4. Values of q s (35 GeV) and a s (44 GeV) 
derived from 0(q s )+NLLA fits to event shape dis- 
tributions, compared with OPAL results. The exper- 
imental and statistical uncertainties are represented 
by the solid, the total error by the dashed error bars. 

Evolving the measurements to y/s = M z o the 
results obtained at 44, 35 and 22 GeV transform 
to a combined value of a s (M z o) = 0.122t°;°og 
which is consistent with the direct measure- 
ment at s/s = M z o by OPAL of a s (M Z o) = 
0.117lo;Qog|||, for the same subset of observables. 

Similarities between the main detector com- 
ponents of JADE j| and OPAL ||, as well as 
between this analysis and studies performed by 
OPAL HjKj at = 91.2, 133, and 161 GeV, al- 
low for a precision test of the energy dependence 
of a s between y/s — 22-161 GeV. The a s values 
from OPAL and from this analysis are shown in 
Figure |. The result of a X 2 fit of the O(o%) QCD 
prediction to the data is shown by the solid line 
corresponding to a s (M z o) = 0. 1207 ± 0.0012 and 
X 2 /d.o.f. = 4.9/5 , taking into account only sta- 
tistical and experimental uncertainties. 
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Figure 5. Values of a s from C(a^)+NLLA fits 
at JADE and OPAL. The solid error bars are the 
statistical and experimental uncertainties, the dotted 
error bars are the total errors. The solid line and the 
shaded band represent a QCD-fit to the data and its 
statistical error, respectively. 

A fit for the hypothesis a s = const, gives a 
X 2 /d.o.f. = 101/5 which has a vanishing proba- 
bility. The energy dependence of a s is therefore 
significantly demonstrated by the results from the 
combined JADE and OPAL data. 
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QUESTIONS 

D. Dissertori, CERN: 

Why did you use only the ln(i?)-matching 
scheme? Typically a set of different matching 
schemes is used, and this in the end could increase 
your systematic error. 

P. A.M. Fernandez: 

In order to compare the present measurements 
consistently with the results given by OPAL, I 
applied only the ln(i?)-matching scheme because 
it is commonly used in OPAL analyses. It has 
turned out that this matching scheme gives the 
best description of data and that the results 
from other matching schemes with comparable 
quality do not significantly differ from the stan- 
dard results. 

L. Trentadue, Parma: 

Is the quality of the data sufficiently good that 
the hadronisation effects are small enough at low 
energies? 

P.A.M Fernandez: 

Of course, the hadronisation effects and its un- 
certainties get larger at lower energies, but they 
remain small enough to perform reasonable mea- 
surements within the fit ranges chosen in this 
analysis, even at y/s = 22 GeV. 

L. Trentadue, Parma: 

It would be useful to consider the variable 
'energy-energy-correlation' for which the analy- 
sis of PETRA data has been done at NLLA level 
in the early eighties. 

P.A.M Fernandez: 

I agree, but I focussed only on the standard 
set of observables which OPAL uses for a s mea- 
surements at high energies, in order to provide 
comparable results over a large range of the e + e _ 
continuum. 



